The nuclear moment of inertia and spin distribution of nuclear levels 
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We introduce a simple model to calculate the nuclear moment of inertia at finite temperature. 
This moment of inertia describes the spin distribution of nuclear levels in the framework of the spin- 
cutoff model. Our model is based on a deformed single-particle Hamiltonian with pairing interaction 
and takes into account fluctuations in the pairing gap. We derive a formula for the moment of inertia 
at finite temperature that generalizes the Belyaev formula for zero temperature. We show that a 
number-parity projection explains the strong odd-even effects observed in shell model Monte Carlo 
studies of the nuclear moment of inertia in the iron region. 
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I. INTRODUCTION 

The shell model Monte Carlo (SMMC) method has 
proven to be quite accurate for calculating the nuclear 
level density in the range of excitation energies up to sev- 
eral tens of MeV [1011]. The advantage of the SMMC 
method is that it can be used to calculate thermal ob- 
servables in model spaces that are orders of magnitude 
larger than can be treated by conventional diagonaliza- 
tion methods. In practice, most of the SMMC calcula- 
tions are carried out in truncated spaces (e.g., one major 
shell) , hence the limitation on the excitation energy. Cor- 
relations become less important at higher temperatures, 
and the results of the truncated SMMC calculations can 
be extended to higher temperatures or excitation ener- 
gies by taking into account the effects of a larger space 
in the independent-particle model 0]. 

The SMMC method can also be used to calculate the 
distribution of nuclear spins at finite temperature 
However, this requires spin projection, and the associ- 
ated computational effort is rather large. For general 
purposes such as constructing tables for large numbers 
of nuclei, simplified models are thus invaluable. The aim 
of this paper is to construct and study such a simple 
model that can reproduce well the spin distributions of 
the microscopic SMMC method. 

A common assumption 0, m global parameteriza- 
tions of nuclear level densities is that the spin distribu- 
tion follows the spin-cutoff model 
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(1) 



Here p is the total level density counting all M states 
in a J multiplet, while pj is the density of spin- J levels 
without the 2 J + 1 degeneracy factor. Thus J2j(2J + 
l)pj = p. The parameter a is known as the spin cutoff 
parameter. The quantities p,pj, and a are all functions 
of excitation energy. 

The model can be derived assuming that the individual 
nucleon spins add up as independent random vectors Q , 



J = 3ii leading to a Gaussian distribution of the spin 
vector P{J) oc e~ J "/ 2cr2 . Integrating over the orientation 
of J, we have P(J) oc J 2 e~ ,/2 / 2cr2 , where J is the mag- 
nitude of the angular momentum (the pre-exponential 
factor J 2 comes from the Jacobian in spherical coordi- 
nates). We recover Eq. by making the semiclassical 
substitution J — ► J + 1/2; the spin cutoff parameter is 
then given by 
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{J 



T2\ 



(2) 



In thermal ensembles it is common to define an effective 
moment of inertia I by the relation between (J 2 ) and 
temperature T, which we can write as 



I 



T 



(3) 



In many of the empirical parameterizations, a is deter- 
mined by this formula using for I the rigid-body moment 
of inertia, I = 2mA(r A 1/3 ) 2 /5, where r « 1.2 - 1.3 
fm is the usual nuclear radius parameter, A is the mass 
number and m is the nucleon mass. Other treatments 
of a based on the independent-particle model have also 
been proposed [9|. SMMC calculations of nuclei in the 
A ~ 50 — 70 mass region show that the assumption of a 
rigid-body moment of inertia breaks down at low excita- 
tion energies starting somewhat below the neutron sepa- 
ration energy, especially in even-even nuclei. The effect 
has a clear odd-even mass dependence. Furthermore, at 
the lowest excitations, deviations are observed from the 
spin-cutoff model itself, and odd-even staggering effects 
(in spin) can be seen. Here we will show that a fairly sim- 
ple model based on a fixed deformation and a fluctuating 
pairing field reproduces very well the detailed SMMC 
results for the effective moment of inertia at finite tem- 
perature. In particular, odd-even effects observed in the 
microscopic SMMC calculations are nicely reproduced by 
a number-parity projection method |10L IllL Il2l Il3 | . We 
would therefore advocate this model for global calcula- 
tions of the spin distributions below the neutron separa- 
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tion energy. Such distributions are needed for theoreti- 
cal estimates of nucleosynthesis reaction rates 7] , among 
other applications. 

Our model is based on the static path approxima- 
tion 0H Gil to the BCS Hamiltonian BCS the- 
ory is valid in the limit when the mean level spacing is 
much smaller than the pairing gap. However, this con- 
dition does not hold in the finite nucleus, in which case 
fluctuations must be taken into account. A similar situa- 
tion occurs in ultra-small metallic particles whose linear 
size is smaller than ~ 3 nm [l^. Theoretical studies 
have indicated that pairing correlations in the crossover 
from BCS to the fluctuation-dominated regime are man- 
ifested through their number-parity dependence. Odd- 
even effects that originate in pairing correlations were 
found in the SMMC heat capacity of nuclei 19]. Such 
effects were also observed in the heat capacities of rare- 
earth nuclei that were extracted from level density mea- 
surements |20L |2jJ. Finite-temperature pairing correla- 
tions at a fixed number of particles were also studied in 
Ref. [13. 

In this paper, we first discuss in Section ITT1 general as- 
pects of calculating the thermal moment of inertia and 
projection on number parity. We work in a grand canon- 
ical ensemble, but the odd-even effects can be extracted 
by the number-parity projection operator. In Section 
IIIII we apply the formalism to a model Hamiltonian that 
includes a deformed single-particle field and a pairing in- 
teraction treated in the static path approximation. This 
yields a formula for the moment of inertia that is a gen- 
eralization of the Belyaev formula [2j| for zero temper- 
ature, explaining the suppression of the inertia at low 
temperature. In Section IIVI we further generalize the 
moment-of-inertia formula to take into account odd-even 
differences, making use of the number-parity projection 
operator. In Section we apply the model to nuclei in 
the iron region using the pfgg/2 shell with single-particle 
energies and wave functions determined from a deformed 
Woods-Saxon potential. The calculated moments of in- 
ertia are found to be in good agreement with the SMMC 
calculations. 



II. FORMAL ASPECTS 



In general, the SMMC method |2J, |25| can be used to 
calculate thermal expectation values of observables O 



(0) = 



Tr(0£ 



-f3H\ 



Z 



where 



Z = Tre 



-0H 



(4) 



(5) 



is the nuclear partition function. H is the nuclear 
Hamiltonian, containing rotational invariant one-body 
and two-body terms. In Ref. 0, 113 exact particle- 
number projection was performed to calculate the traces 
in Eq. Q at fixed number of protons and neutrons. 



In Ref. the spin distribution was calculated us- 
ing spin projection techniques. For temperatures that 
are not too low, it was found that the spin-cutoff model 
(|T|l describes rather well the spin distribution but with 
an energy-dependent moment of inertia. The purpose of 
the present work is to understand the temperature de- 
pendence of the moment of inertia in terms of a simple 
model. We note, however, that at the lowest tempera- 
tures the SMMC calculations reveal deviations from the 
spin-cutoff model JTJ , which are beyond the scope of the 
model discussed here. 



A. Moment of inertia 

In this work we shall assume that the spin distribution 
can be described by Eq. , and we therefore only need to 
calculate the variance a 2 . The obvious way to do this is 
to evaluate the expectation value of the operator O = J 2 
directly from Q}, as is done in SMMC. However, our 
model in Section ITTT1 is based on a deformed Hamiltonian 
-ffdef, and for such Hamiltonians it is useful to define a 
moment of inertia tensor as the response of the nucleus 
to a rotational field Co. 

We shall work in the grand-canonical ensemble, replac- 
ing H by H' = H dci - fiN in Eqs. © and (QJ. In the 
presence of a rotational field, the Hamiltonian is given 
by H' — lu ■ J and its free energy is 



F{(3,uj) = -/rMnTre-^'-^ 



(6) 



The moment of inertia Iij is defined by the expansion of 
F to second order in (3, F(/3, w) = F((3, uj — 0)— 
where u>i are the components of lu. Equivalently 26] 



y dujiduj 



d 2 F f 13 



where 



rrur p -f3H' < e rB' j -tH'\ t.i 

<Ji(r)4(0)> = ^ ^ (8) 

is the spin response function in imaginary time. For 
a rotationally invariant Hamiltonian, Jj(r) = J^, and 
Iij = ISij with I = f3(J 2 ) = /3(oP)/3, in agreement with 
Eqs. and ©. Choosing the cranking axis along the 
fixed z-axis of the laboratory frame, we can calculate / 
from 



- duj2 \u=0 , 



where 



F(f3,uj) = -TlnTre-^ H '- uL ^ 



(9) 



(10) 



1 Here we use a different notation for the angular momentum in the 
laboratory frame to distinguish it from the angular momentum 
in the intrinsic frame. 
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and the angular momentum component along z is de- 
noted by L z . 

A non-rotational invariant effective Hamiltonian arises 
in the mean-field approximation when the single-particle 
potential is deformed. In such a case Pdef describes the 
Hamiltonian in the intrinsic frame of the nucleus. The 
quantity /y in Eq. is then the moment of inertia 
tensor in this intrinsic frame, where Ji are the intrin- 
sic components of the angular momentum J. To recover 
the moment of inertia I in , it is necessary to integrate 
over all orientations of the intrinsic frame and then use 
@. One obtains the result (see Appendix B) 



1 , 



g yuiX * *yy 



hz) ■ 



(11) 



Eq. expresses the effective moment of inertia / in 
terms of the intrinsic principal moments In. 



where we have used the notation 

Tr (Oe^e-P 11 ' 



Tr 



N p -0H> 



(16) 



The number-parity projected moment of inertia R is 
defined from the second-order expansion (in lu) of the 
number-parity projected free energy 



F„ 



-TlnTr 



We find 



Xf dr {Mr)J 3 {Q))+v{e" N )^dr (J,(r) J,(0)). 



p e -f3{H'-Z-J) 



(17) 



TV - 



1 + r)(eJ 



where {Ji{r)Jj{Q))^ is defined as in (|To|) . 



(18) 



Number-parity projection 



III. MODEL 



The calculations in the previous Section III Al were de- 
scribed in the grand-canonical ensemble. While this al- 
lows the average number of particles to be specified, it 
is not precise enough to reproduce odd-even effects. We 
note that the behavior of odd and even nuclei at low 
temperatures is quite different; the spin goes to zero 
for even nuclei due to pairing, but remains finite at 
zero temperature for odd nuclei. Exact particle-number 
projection can be done using the projection operator 
P N = e l ^ N -^d(f>/2Tr as in the SMMC, but leads to 
cumbersome expressions. In order to capture the main 
odd-even effects, it is often sufficient to use a number- 
parity projection |lCt llll ll2L 1 1 3| that distinguishes only 
between even and odd number of particles. The number- 
parity projection is defined by 



rje 



(12) 



where r\ = 1 or -1 describes the projection on an even or 
odd number of particles, respectively Thus, the number- 
parity projected partition function is 



Tr (P n e"^') 



1 



(13) 



where the bracket denotes a thermal trace, (O) = 
Tr (^Oe~P H /Tre^^H . We can also calculate number- 
parity projected expectation values of observables 



(0)r,= 
Using ljT2|) . we find 



Tr ( OP n e~ m ' 



Tr {Pr,e-P H ') 



(14) 



We now ask, starting from the independent-particle 
shell model, what is the minimal model that will include 
the most relevant interaction effects for calculating the 
spin distribution. Clearly, the most important correla- 
tions are those associated with the quadrupole deforma- 
tion and the pairing field. Both of these can be treated 
in a mean-field approximation, but the mean-field equa- 
tions predict sharp transitions that are not supported by 
more detailed theories. Thus we go one step further in 
the finite-temperature theory, using the static path ap- 
proximation (SPA) Q m 1 1 6| of the partition function 
to include time-independent fluctuations of the order pa- 
rameters. 

We consider an Hamiltonian composed of an axially 
deformed Woods-Saxon well for the single-particle poten- 
tial and orbital-independent pairing for the interaction. 
We denote by | k) the single-particle eigenstates in the de- 
formed potential with energies ek- They can be divided 
into degenerate time-reversed pairs (k, k). For an axially 
symmetric potential, \k) — \q, /i) where [i is the pro- 
jection of the angular momentum on the symmetry axis 
and q are other labels of the states. The time-reversed 
states are defined by \k) — \q, — /i) = | — k) (known as 
the BCS phase convention), and we adopt the convention 
k > /i > 0. The Hamiltonian may then be expressed 
in the form 

ffdef = ^{a\a k + 4 as ) - GP^P , (19) 

fc>0 

where is the pair creation operator, pt = J2k>o a l a I> 
and G is the pairing strength. 



A. Static path approximation 



The Hamiltonian (|19|) contains a pairing interac- 
tion. Using the Hubbard-Stratonovich transformation, 
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the imaginary-time propagator e~^ H can be written as 
a functional integral over pairing fields of propagators 
that describe non-interacting quasi-particles. Here we 
shall use the SPA, which takes into account only time- 
independent pairing fields. The functional integral then 
reduces to an ordinary integral over a complex pairing 
field i 



-0H> 



(3G_ 
2^ 



(20) 



where 



H k = (e k -n-G/2)(ala k +ala- k )-GCa- k a k -G^alal+G/2 

(21) 

Our model (|19(l describes nucleons moving in a de- 
formed well, but it could have been derived from a rota- 
tionally invariant Hamiltonian that included quadrupolar 
two-body interaction. This would introduce five addi- 
tional integration variables in the SPA integral (|2UI) , two 
representing the intrinsic deformation and three repre- 
senting the orientation of the deformed field 0| . This 
integration over the Euler angles of the intrinsic frame 
is equivalent to the symmetry restoration described by 



1. Partition function 

Using H20|) , we can represent the grand-canonical par- 
tition function in the form 



Z ^ I 

2n 



Here we used Tr (J\ k 




fee 



-0H h 



(22) 



fe>0 



Ti-fc is the trace evaluated in the Fock space of the or- 
bital pair (fc,fc), i.e. in the 4-dimensional space spanned 
by {\n k: n k )} = {|0, 0), |0, 1), |1, 0), |1, 1)}. In this repre- 
sentation, H k is the matrix 



H k = 



(G/2 



\GC 





e k - n 

e k - fi 







2(e fe - M )-G/2/ 



(23) 



The traces in l|22|) are easily evaluated by diagonalizing 
each H k in the corresponding 4-dimensional space. The 
four eigenvalues are e k — [i + {~E k , 0, 0, E k }, where 



E k = v/(efe-M-G/2) 2 + G 2 |e| ; 



(24) 



are the familiar quasiparticle energies 2 but now defined 
for an arbitrary complex pairing field £. The trace in the 



2 One usually denotes the self-consistent G|§| as A in BCS formu- 
lation. 



subspace (fc, k) is then easily evaluated as 

Tr k e~ l3Bk = e- l3 ^"- fi Hcosh 2 ((3E k /2) 

— e -P(ek-v)(i _i_ e ~' 3£; fc)(l 4- e l3Ek ) 



(25) 



The last algebraic form is convenient when dividing by Z 
in the evaluation of expectation values, as the reciprocal 
is proportional to f k (l—f k ) where f k are the quasiparticle 
occupation probabilities. 

An alternative way of calculating the trace is to write 
H k = (e k -n) + (4 a- k )H k (^) where H k is the 2 x 2 
matrix 



and use the identity 0] 



efe M 2 "s ] (26 ) 



Tr exp 



(4 a k )JC i a J 



det (l + e K ) 



(27) 



for the matrix K. = —fiTL k . The eigenvalues of 7i k are 
just ±E k , leading again to Eq. |2HJ)- 

The complete grand-canonical partition function is 
given by 



Z = B. I dtdee-™ -^^-") TT 4cosh 2 (/3£: fe /2) 
2n J fe>0 

(28) 



It can also be written in the form 



2n 



(29) 



Tr k e-^ H ", where where 



F K, C] = G|e| 2 +^(e fc -^)-2/3- 1 E ln I 2 cosh(/3S fc /2)] 

fc>0 fc>0 

(30) 

is the free energy for a complex pairing field £. 

The canonical partition function Zjy = Trjye - ^ 
can be calculated by a Fourier transform of the grand- 
canonical partition 



(31) 



Z N = — / da e~ aNr Tr{e aN e~ fm ) 
2m J_ l7T 



The integrals in (|31|) can be evaluated in the saddle 
point approximation in both £ and /i. The corresponding 
saddle-point equations, dF/d^* = and dF/dfi = —N, 
will give the usual finite-temperature BCS equations 



G Z - 

fe>0 



tanh 



2E k 



(32) 
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and 



N = 



E 

k>0 



1 - 



£fe - /X- 



■ tanh 



(33) 



where the quasi-particle energies -E/c are given by (|24H . 
The solutions of (|32l) and (|33[) determine the pairing gap 
A = G|£| and the chemical potential /i as a function of T 
and particle number ./V (the phase of £ is undetermined) . 

A better estimate of the canonical partition function 
Zn can be obtained by a saddle-point integration in i|31|) 
over a for every £, but keeping the integration over £ 
intact. We find 



Z N « 2/3G / d|C| |f| 2^T 



ctyi 2 



-1/2 



-/3(^K,r]+^) 



(34) 

where n = /i(iV, T, £) is determined from dF/d/j, — —N, 
i.e., Eq. (J33J, and 



d 2 F 



^ /?£ fc (e fc - fj, - ±Y + sinhQ3£ fc )G 2 |£h 

fe>0 



2E 3 k cosh 2 



(35) 

In Eq. (|34ll we have carried out explicitly the integral 
over the phase of the pairing field £ since the integrand 
was only a function of |£|. 



2. Moment of inertia 

Our model (|19|) describes a non-rotationally invariant 
Hamiltonian, and we can use the formalism of Section 
111 Al to estimate the moment of inertia I in terms of the 
intrinsic moments la (see Eq. Rather then using 

Eqs. and JHJ with the full Hamiltonian H, we first 
apply an SPA representation similar to (|20(l but for the 
cranked Hamiltonian in 10, and then calculate the in- 
trinsic moments from la = d 2 F(P,uj)/doj 2 \ LU= o- If this 
is done starting from the canonical partition function in 
(|6*|) , we obtain the following expression 



J d£d£* _i_ J*^ da e- aN e-mw] 
where is the free energy l|3Ufl and 



(36) 



2y(£)= /■'dr'&[e-' Jir «J i (r)J J -(0)]/'IV[e-^], (37) 
Jo 

where Ji(r) = e TH ^Jie~ TH i and = Y,k>o H k is tne 
mean- field Hamiltonian in a pairing field £. Expression 



(|37|l is analogous to Q and JSJ), except that the Hamil- 
tonian H' in those expressions is now replaced by H^. 

The integrals over a in i|3fi|l can be done in the saddle- 
point as before to obtain the final expression for the in- 
trinsic moments 



Iij — 



S" m |g|(a 2 F/V)- 1 / 2 e-^K-e , ]+^)z ij .(g) 
/ °° d \Z\ \£\(B 2 F/dn 3 )- 1 / 2 e-M F foW+* N ) 



It remains to calculate the moments Iij (£). The opera- 
tor J z leaves the (fc, k) subspace invariant, but the opera- 
tors J x and J y connect different subspaces k and k' , so the 
trace in l|37|) is to be evaluated in a 16-dimensional space. 
This is most conveniently done in the quasiparticle rep- 
resentation. The transformation from deformed single- 
particle states to the quasi-particle states is achieved 
through a Bogoliubov transformation 



a J)=u k ( a j 



Uk —Vk \ I a k 



(39) 



where is real and Vk is complex, and u\ + \ v^\ 2 = 1 to 
preserve the fermionic commutation relations. Relations 
(HSU imply 



= Uk ; fj, = — Vk for k > 



(40) 



The parameters Uk , Vk are chosen such that Hk in 
Eq. (|21|l is diagonal in the quasi-particle representation, 
i.e., 



(41) 



where 7ifc is the 2x2 matrix (|26[) and Ej. are the eigen- 
values 124|) of 7ife . The solution is 





e/s - /x - 


G/2 






K 1 - 


£fe - M - 


G/2 







(42) 



and argufc = arg£. The Hamiltonian iij, is now given by 
#fc = (efc - A*) + £fc(c4 a fc - ot^.a\), and 

H e = £ E ^ a W + a i a k) + E( £fc -v-Ek). (43) 



fc>0 



fc>0 



Expressing Ji in the quasi-particle representation, and 
using Q43|l ■ we can calculate the intrinsic moments in 
closed form (see Appendix A). The final result is 
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fc,;>o 

(44) 

where 

are the quasi-particle occupations, Uk and v k are still given by (|42|) . except that now we have chosen 

u k ,v k > for k > , (46) 

and u%, v k are still given by 1)41) [I. Eq. 144|) is the finite-temperature generalization of the Belyaev formula |23| . 

Eq. I|44|) can be rewritten by separating out the contribution from the k — I terms in the sum. Using (/; — fk)/ (E k — 

Ei) -dfk/dE k , we obtain (for i = j) 

i u (o = 2j2(\m\k)\ 2 + \m\ fc)n + 2 E (imoi 2 + \m\ - oi 2 ) 

fc>o ^ k ' k,i>o (47) 

^ f / | ^2 fl - fk , , n2 1 _ /* - // ^ 
X <^ (lifeM; + Ufe^i) — — + (UkVl -VkUl) — —— 

{ &k — rj\ &k + &l 

In particular, (k\j z \l) — for k =/= I, and (k\j x \k) — 0. Therefore, the moment of inertia around an axis parallel to the 
symmetry axis z (non-collective rotation) is given by 

X zz (0=2El<^|fc>| 2 (-f! i ) , (48) 

k>0 ^ k ' 

while the moment of inertia around an axis x perpendicular to the symmetry axis (collective rotation) is given by 

i xx (o= 2EK fc i^i- fc )i 2 (-|^)+ 2 E imm\ 2 +\{m-i)\ 2 ) 

k>a ^ k ^ k,l>o (49) 

„ J / i \2 fl ~ fk . , x2 1 — fk — fl 

X < (UkUi + V k Vl) — — + [U k Vl - V k Ul) — — — — 

^ £jk~ Ej\ £jk + &\ 

In the limit T -> (but A > 0), -df k /dE k -> S(E k ), and since E k > 0, = 0. Also /* ->• and Eq. gSJ 

reduces to 



4,(0 =2 x: owjoi a +iw.i-oi a ) ^"^ =Eiw-io r" fc g" r • (50) 



Eq. H5()|l is known as the Belyaev formula |2jj ; it produces 
a moment of inertia that is suppressed relative to the 
rigid-body value. 



IV. NUMBER-PARITY PROJECTION 

In the relations we derived in the previous Section for 
the partition function and moment of inertia, the number 
of particles is fixed only on average, and odd-even effects 
cannot be reproduced. Here we go through the same 
derivation steps but include now the number-parity pro- 
jection operator P^. The resulting formulas will exhibit 
explicit terms depending on the number parity. 



Partition function 



The projected partition function (|13|) introduces the 



operator e 



ittN 



In the SPA 



Tr U^e-P H ') 



(3G 
27 



-/3G|5| 2 



fc>0 



e M<4 <**+<»£»£) e -/3H k 



• (51) 



Within each subspace (k,k), the operator e t7rN changes 
the sign of the two vectors |0,1),|1,0), but leaves the 
sign of 1 0, 0} , |1, 1) unchanged. The matrix representing 
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o v fc *-r fe a.; lg then 



/l 














-1 














-1 





Vo 








1 



(52) 



The transformation that diagonalizes Hi- in Eq. ll'ol) 
leaves this matrix invariant, and therefore the trace is 
now given by 



Tri. 



e iir(ala k +ata- k ) e -f3H k 



= e -/3(e*,-M)( 1 _ e -0E k ^ e 0E k _ x y (53) 



The projected grand-canonical partition function Z v 
Tr \P n e~P H j is now calculated from and to 



be 



/3G 



1+77 n tanh; 



fc>0 



(54) 

Notice that the integrand in (|54f) has the form of Eq. 
when applied to the Hamiltonian — X)fe>o ^ k a ^ a 
fixed pairing field £. Indeed 



Tr (e^e-ext 



Tr e"^ 



n tanii: 



(55) 



fe>0 



This projected partition can also be written as Z v 
^Jd^dC e -W«>n where 



'l + '/n fc >otanh 2 ^~ 



(56) 

is the number-parity projected free energy. Proceeding 
as m Section IIII A II we can derive (in the saddle point 
approximation) number-parity projected BCS equations 




C 



E k smh(pE k ) 



(57) 



~ (tn,:):-. 



2C n 



smh(/3E k ) 



^IWanh 2 (g 
l + r?n fc>0 taah 2 (^ 



(58) 



(59) 



It is interesting to take the T — > limit for the above 
equations. For ry = +1, they simply become the usual 
T = BCS equations. However, for the odd projection 



rj = —1, we find (assuming there are no degeneracies and 
A = > 0) 



-=r— 



1 



G — ' 2i?fe 1Eh t , 

fc>0 
|2 



fc#fen>0 



1 



(60) 



where |«fc| 2 is given by (|4*2|) and i?fe is the lowest quasi- 
particle energy (corresponding to Ck closest to /i). In de- 



riving l)6()|l. we have used the limit C_i 



/3- 



o/4. 



Eqs. (|60|) reproduce what is known as the blocking ef- 
fect, since one level k = kg is "blocked" and does not 
contribute to the sum over k. 

In Fig. ^ we display the solution of the number-parity 
projected BCS equations for a Hamiltonian correspond- 
ing to the nucleus 56 Fe. The pairing gap A is shown as 
the solid line as a function of temperature T . For the 
particle- number projected BCS equations I|57J1 . the solu- 
tions for even and odd particle numbers are shown as the 
dot-dashed and dashed lines, respectively. The proton 
gap A p and the neutron gap A„ are shown in the left 
and right panels. Note the strong suppression of the gap 
for the odd projection. Our results for the projected gap 
are qualitatively similar to Ref. |12lll3 |. but the formulas 
are quite different. 
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FIG. 1: The pairing gap A versus temperature T for protons 
(left) and neutrons (right). The solid lines are the solution of 
the BCS equations 132H for 56 Fe. The dotted-dashed lines and 
the dashed lines are respectively the solution to the number- 
parity projected BCS equations 1571 for rj = 1 and rj = — 1. 
The deformation parameter is taken to be /?2 = 0.24 and the 
pairing strengths are determined from the zero-temperature 
BCS to reproduce the experimental values of the gaps (using 
the second difference formula of Ref. 12711. 



The number-parity projected partition function at a 
fixed average number of particles is given by an equation 
similar to (|34f> except that F is replaced by F v . 
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Moment of inertia 



Thus we need to calculate the projected value 
( Jj (t) Jj-(0))„. The odd-even number projection can be 
The number-parity projected moment of inertia can be carried out for any operator O using (jEU) at a fixed pair- 



calculated as in Section llTTX2l but now starting from the in S field £ together with ^ 
number-parity projected free energy in the presence of a 
rotational field uj. The result is 



13 f™d\£\ \£\(d 2 F ri /dv?)- 1 / 2 e-P( F n[^*]+^) ' 

(61) 

where 



(0) v 



l + ^(lL>o^nh 2 ^) 



(63) 



w= 



dr 



Tr [P v e-f )H tJ i (T)J j (0)] 
Tr [P v e~P B i] 



In general, (O)^ = Tr IOe inN e~^ j /Tr (e iirJV e -/JH « 

will have the same form as (O) but with e ±/3Ek replaced 
(62) by-e ±(3E *. 



For the intrinsic moment of inertia we find 

dr (JitfJj (0)) + ?7 n fc>0 tanh 2 fiffc dr (Mr) J, (0)% 



l + ^n fc>0 tanh 2 ^ 



(64) 



where J Q dr (Ji(r) Jj(0)) is given by l@3J and J P dr (Ji(r)J J (0)) 7r is obtained from the expression for dr (Jj(r)Jj(O)) 
by the substitution 



/fc — > /fc — 



1 



1 _ e 0E k 



(65) 



We now inspect the T — > limit of Eq. (|64|l . For the even number-parity projected inertia ZJ^ -1 — ► and Z^r 1 is 
the same as in Eq. I|5UI) . For the odd number-parity projected moment of inertia, we have in the limit T — > 



2'J 



/3|Mo| 



(66) 



where |fco) = \q, fJ-o) {Ek is the lowest quasi-particle energy), and 
l^' 1 - 2 £ (| < fc|j x |fc > I 2 + I < %*| - fc > | 2 ) 

fc^fc o >0 

+2 ^ (| < k\j x \l > | 2 + | < k\j s \ - I > | 2 

k,l^k >0 



(ukUkp + VkVko) 2 (u k Vk - VkUkpf 

E k - E k( 
(u k vi - V k Ul) 2 



E k + Ei 



E k + E ka 
+ P\ < fcob'xl -ho > I 2 



(67) 



The number-parity projected moments of inertia 1^ 
are computed by a numerical integration of Eq. I|61|l , us- 
ing Eq. H44fl and (|64|l . For the final result, the contri- 
butions from the three principal axes must be combined 
according to Eq. (JTJJ. 



V. COMPARISON TO SMMC RESULTS 

We have used our formulas to study the moment of 
inertia I versus temperature for nuclei in the iron region 
and compared the results with SMMC calculations. The 
effective Hamiltonian is defined in the full pfgg/2 shell 
with single-particle energies determined by a Woods- 
Saxon potential plus spin-orbit term. The interaction in- 
cludes T = 1 monopole pairing and multipole-multipole 



interaction terms (up to hexadecupole) pj . In SMMC we 
calculate (J 2 ) as an observable and then find I from 

is = P( J *) ■ ( 68 ) 
n 

Significant odd-even effects are observed in the SMMC 
calculations see, e.g., in Fig. [21 where the SMMC 
thermal moment of inertia of the even-even nucleus 56 Fe 
(solid circles) is compared with the SMMC moment of 
inertia of the odd-even nucleus 55 Fe (open circles). We 
see stronger suppression at the lower temperatures in the 
even-even case. 

The parameters of the model we discuss in the present 
work are the deformation fa and the pairing strengths 
G p , G n for protons and neutrons, respectively. We have 
chosen fa — 0.14 and computed G p and G n from the ex- 
perimental values of A p and A„ (using zero-temperature 
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FIG. 2: Moment of inertia (characterizing the spin-cutoff 
distribution versus inverse temperature /3 for an even- 
even nucleus ( 56 Fe) and an odd-even nucleus ( 55 Fe). The 
SMMC results calculated from I68H are shown by the solid 
circles for 56 Fe and the open circles for 55 Fe. The results of 
the model discussed in this work (which includes fluctuations 
of the pairing field and number-parity projection) are shown 
by solid line for 56 Fe and dashed line for 55 Fe. For protons 
we use even number-parity projection, while for neutrons we 
use even (odd) number-parity projection for 56 Fe ( 55 Fe). In 
the model we use a deformation of @2 = 0.14 and pairing 
strengths of G p = 0.42 and G„ = 0.36. 



BCS). In our model, we include fluctuations in the pair- 
ing gap and we have found it necessary to reduce the 
BCS values of G by about 20% to reproduce correctly the 
high temperature behavior oil . We then fixed /?2 = 0.14, 
G p = 0.42 and G n — 0.36 through the whole mass region, 
and calculated the number-parity projected moment of 
inertia T" 1 (with fluctuations in the pairing order param- 
eter included) using Eqs. Jill), ijrTT)!. and (jSJ. Fig.0 
shows the results for I v in 56 Fe (solid line) and in 55 Fe 
(dashed line). We find that the inclusion of a number 
parity projection allows us describe reasonably well the 
odd-even effect observed in the SMMC moment of iner- 
tia. We note that an exact particle-number projection 
(for both protons and neutrons) is used in the SMMC 
method. 

To demonstrate the importance of fluctuations in the 
pairing order parameter, we show in Fig. El the number- 
parity projected moments of inertia I n (£) evaluated at 
the corresponding solutions £ to the number-parity pro- 
jected BCS equations l|57(l but without the inclusion of 
fluctuations. I v {£,) is calculated from l|rT4"| and averag- 
ing over principal axes. We use the same deformation 
and pairing parameters as in the theory with fluctua- 
tions. These results demonstrate that we cannot obtain 
good agreement with the microscopic SMMC calculations 
without including fluctuations. 

In Fig. 01 we present a systematic study of the moment 
of inertia for even and odd iron isotopes from 55 Fe to 



FIG. 3: Similar to Fig. but the lines correspond to the 
number-parity projected moments of inertia calculated at the 
solutions £ of the corresponding number-parity BCS equations 
1571 . The deformation and pairing parameters are the same 
as in Fig. [5] 

60 Fe. We show the moment of inertia versus /?, com- 
paring the SMMC results (symbols with statistical error 
bars) to our model results (solid lines for even isotopes 
and dashed lines for odd isotopes). The model includes 
fluctuations in the pairing fields, and we have used a fixed 
set of deformation and pairing parameters. Despite its 
simplicity, the results of our model agree well with the 
full microscopic SMMC calculation. 
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FIG. 4: Systematics of the moment of inertia versus f3 for 
a series of iron isotopes. The symbols (solid circles for even- 
mass isotopes and open circles for odd-mass isotopes) are the 
SMMC results calculated from 1681 . The lines describe the re- 
sults of our model (including fluctuations in the pairing field). 
The solid lines are the even number-parity projection while 
the dashed lines describe the odd number-parity projection. 
The dotted-dashed lines are the rigid-body moment of inertia. 
We use /3 a = 0.14 and G p = 0.42, G n = 0.36. 
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FIG. 5: As in Fig.0] except that the moments of inertia are 
shown versus temperature. 
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FIG. 6: The moments of inertia versus temperature for the 
odd-even nucleus 55 Mn and the odd-odd nucleus 56 Mn. The 
notation as in Fig. [I] 

In Fig. [5] we show the same results as in Fig.0]but plot- 
ted as a function of temperature. In general, the thermal 
moment of inertia is seen to decrease with decreasing 
temperature with stronger suppression in the even iso- 
topes. Furthermore, for the odd iron isotopes (in which 
an odd number-parity projection is used for neutrons), 
the moment of inertia starts to rise at low temperatures. 
Indeed, according to (|fj5|l 

I -> ~/3J(J + 1) , (69) 

for large (3 where J is the ground-state spin. In an odd- 
even nucleus, J ^ and / increases linearly with f3 at 
large (3. In the model this rise at low temperatures can 
be understood to be the effect of the unpaired neutron. 
Using the projected moment of inertia for odd number of 
particles, we find in the limit of large (3 

/^-^(Mo + l^olixl-fco)! 2 ), (70) 



where the second term in (|70|) contributes only for /j,q = 
1/2. It is interesting to note that a similar odd-even effect 
was found in the spin susceptibility of ultra-small metal- 
lic particles (nanoparticles) [28|. The rise of the spin 
susceptibility at low temperatures for an odd number of 
electrons is known as re-entrant behavior and survives in 
the fluctuation-dominated regime. 

Finally, in Fig. we show the moment of inertia versus 
temperature for an odd-even nucleus 55 Mn (left panel) 
and for an odd-odd nucleus 56 Mn (right panel). We have 
used the same deformation and pairing strength parame- 
ters as for the iron nuclei. Again we find good agreement 
between the results of the model (dashed lines) and the 
microscopic SMMC calculations (symbols). 



VI. CONCLUSION 

We have presented a simple model to calculate the nu- 
clear moment of inertia at finite temperature. The model 
includes quadrupolar deformation of the single-particle 
field and pairing interaction. The pairing interaction is 
treated beyond the mean-field BCS limit by including 
static fluctuations in the pairing order parameter at fi- 
nite temperature. 

In the fluctuation-dominated regime, finite tempera- 
ture signatures of pairing correlations are often observed 
as odd-even effects (in particle number) 0. Such sig- 
natures are usually not seen in the grand-canonical en- 
semble, but only after particle-number projection is im- 
plemented (such as in SMMC). In this work we have 
used number-parity projection instead of exact particle- 
number projection and showed that the odd-even effects 
in the moment of inertia seen in the microscopic SMMC 
calculations can be well reproduced. The advantage of 
such number-parity projection is its simplicity, which al- 
lows us to express the projected moment of inertia in 
essentially closed form (except for an integral over the 
pairing field). 

The simple model developed here is useful in estimat- 
ing the spin distribution of nuclear level densities at not- 
too-low temperatures for which the spin-cutoff model 
usually works well. The spin-cutoff model depends on 
a single parameter - the moment of inertia, and thus 
our model is useful for global estimates of the spin dis- 
tribution below the neutron separation energy where the 
moment of inertia may deviate from its rigid-body value. 
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Appendix A: finite-temperature generalization of the Belyaev formula for the moment of inertia tensor 



Here we derive Eq. H44jl. Since is diagonal in the quasi-particle representation (see Eq. I|43[) ). we have 

4(t) = e TH ia>le- TH t = e TEk a{ 
a k (r) ee e TH ta k e- TH z = e- rEk a k . 

Using (|71|) and the Bogoliubov transformation l|39|) . we obtain 
Mr) = Um\lW k (r)ai(r) =^(k\j i \l)(u k ai(T)+v k a- k (T))(u l a l (r)+v l aUT)) 

k,l k,l 



k,l 



'u k uia' k ai 

The spin response function is then given by 
(Ji(T)J,(0)) = £ miQimljjln) x r 

+ e r(iS fc -£; 1 ) UfeW;UmWn ( a t a/Q , j _ iQ ,tj + e -r(B fc - J B i ) 



(71) 



(72) 



A'.?. 



Ufc "m v n (a k a\ am a' n ) 



+ e T{ - Ek+El) u k viv m u n {a\a\a 7 na n ) + e T{Ek+El) v k uiu m v n (a k aia m a\) 
Using Wick's theorem in the quasi-particle representation and f k = f k , we find 

(Mt)J 3 (0)) = £){<fc|*|fc>[t£/ fc + v 2 k (l - f k )]} x J2i( m \3oH[u 2 m fm + v 2 m (l - f m )}} 



(73) 



J2(k\Ml)(l\m X [ e ^ £ - £ ') U 2 fe ^ /fe(1 _ fl) + e -r { E k -E l)v 2 v 2 {l _ fk)fi 

+ e< E *+ E Mv?hfl + e-^ E ^vluf{l - f k )(l ft) 



k.l 



x u k v k uivi 



r(E k -Ei 



- f,\ - p - r (Ek-E l ) ( , _ 



A(l-/0-e 



(1 - fk)fi 



+ e^+^hh + e-r{E k +E l){1 _ /fe)(1 _ fi) 



(74) 



(75) 



The sum in the first term is zero by symmetry (we assume axially symmetric deformed states). To proceed, we need 
to carefully examine the behavior of the matrix elements of ji under time reversal. Denoting by T the time reversal 
operator (cf. Ref. j^, l-2c), we have TjiT^ 1 = — jj, and for k > (i.e., for spin projection [i > 0) 



T\k)=T\q, ri = \q, —fi) = \k) 
T\-k)=T\q, -/i) = -|g, /x) = -|fc> . 

Using the transformation properties of the matrix elements under time reversal [see (1-34) in pTj]. we then have 



(76) 



{k\ji\l) = -(l\ji\k) for k-l>0 

mm = mm for k-i<o. 

Using relations 177|) . we can rewrite l|74|l in the form 

(j<(t)4(q)) = e imm\jj\k) + (i\m(k\ j3 \i) + mm - 0Hii#) + - *» 

fc,Z>0 

x [e T ( B *- B »)( tt 2 tt 2 + UkVkUm )f k {\ - /,) + e -^ Ek - E ^(v 2 k vf + u k v k uivi)(l - f k )f t 

+e T(E k+El )^ u 2 kV 2 _ UkVkUlVl )J k fl + e -r(B fc +B,) (w 2 u 2 _ UfcWfcUlUl )(l _ / fe )(1 _ /,)] 

The intrinsic moments are given by 2y = J Q dr (Ji(r) Jj(0)). Using (|78|) and the relations 



(77) 



(78) 

(79) 
(80) 



f k (l - fi) f e< E ^dr = (1 - f k )fi f e-^- E 'Ur = fl fk 
Jo Jo 

fkfi f ' e- T ^ +E ^dr = (l-f k )(l-f l ) f \r{E k+ E l)dT = 
Jo Jo 



E k -Ei 

1 ~ fk — fl 

E k + Ei 



(81) 



we arrive at Eq. (|44|) . 
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Appendix B: derivation of the effective moment of inertia 

Here we derive Eq. (|llf> . The components of the angular momentum in the laboratory frame are related 
to the components Ji in the intrinsic frame through = / 2?*^], (ip, 9, (fyJ^' , where T>^ are Wigner rotation 
matrices and (ij),0,<f>) are Euler angles. In particular L z — '^2, i u)iJi where the unit vector ui is given by u> = 
(— sin^cos^, sin#sin</>, cos(9). To restore rotational invariance, we integrate over all the possible orientations of 
the intrinsic frame. In practice we replace in (|10J) 



r d ^ f 09 sm6e-^ H '- u ^ aiJi \ (82) 
4tt Jo Jo 



e 



and use ©to find 



I = — J d(f> j d9 sin 9'^^ujilijilij . (83) 



where ly is the intrinsic moment of inertia tensor (JHJ). The explicit integration of Ij83(l leads to (|11J) . 
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